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Abstract 

The fundamental problem of the transition from quan- 
tum to classical physics is usually explained by decoher- 
ence, and viewed as a gradual process. The study of en- 
tanglement, or quantum correlations, in noisy quantum 
computers implies that in some cases the transition 
from quantum to classical is actually a phase transi- 
tion. We define the notion of entanglement length in d- 
dimensional noisy quantum computers, and show that 
a phase transition in entanglement occurs at a critical 
noise rate, where the entanglement length transforms 
from infinite to finite. Above the critical noise rate, 
macroscopic classical behavior is expected, whereas be- 
low the critical noise rate, subsystems which are macro- 
scopically distant one from another can be entangled. 

The macroscopic classical behavior in the super- 
critical phase is shown to hold not only for quantum 
computers, but for any quantum system composed of 
macroscopically many finite state particles, with local 
interactions and local decoherence, subjected to some 
additional conditions. This phenomenon provides a 
possible explanation to the emergence of classical be- 
havior in such systems. A simple formula for an upper 
bound on the entanglement length of any such system 
in the super-critical phase is given, which can be tested 
experimentally. 



1 Introduction 

Quantum computation is a fascinating subject which 
manifests the peculiarities of quantum mechanics and 
uses them in order to achieve an advantage in terms of 
computational power over classical computers. Shor's 
algorithm ||23| is the most astonishing known example 
for such an advantage in computational power: It en- 
ables one to factor an integer in polynomial time using 
a quantum computer, whereas the best classical algo- 



* Computer Science 
doria@cs.berkeley.edu 



division, UC Berkeley, E-mail: 



rithm for this task is sub-exponential. The advantage 
of quantum algorithms such as Shor's algorithm over 
classical algorithms suggests, though not proves, that 
the computational complexity class of quantum com- 
putation is not polynomially equivalent to that of clas- 
sical, even randomized, Turing machines. 

In real life, quantum computers will be subjected 
to noise. We will make here the assumption that the 
noise is probabilistic and local, meaning that each par- 
ticle, each time step, suffers a certain faulty event with 
probability ?/, which is referred to as the noise rate. 
Quantum computing is now known to maintain its full 
computational power even in the presence of such local 
noise, as long as the noise rate is weaker than a cer- 
tain threshold tjq^, ^, ^. On the other hand, it is 
knownH that when the noise in the system is stronger 
than a much higher threshold, rji, the quantum compu- 
tation can be simulated efficiently by a classical Turing 
machine. Trying to put together the two results, we 
learn that there are two regimes of noise in the quan- 
tum computer, in which the computational power of 
the system is qualitatively different. For weak noise, 
the computational power is fully quantum, whereas for 
strong noise it becomes classical. This raises the fol- 
lowing question: What is the physical difference be- 
tween the two noise regimes, which reflects itself in the 
difference in computational power, and how does the 
transition between the two different physical behaviors 
occur? 

It turns out that an answer to these questions can 
be given in terms of the behavior of entanglement, or 
quantum correlations in the system. Perhaps the best 
way to explain the notion of entanglement is by saying 
what entanglement is not: We say a state in the Hilbcrt 
space of a composite system A(^ B is non entangled, if 
two persons, Alice and Bob, each having access to one 
side of the system, could construct the overall state by 
applying local quantum operations on their side, and 
exchanging classical information, by, say, talking on the 
phone. Any state in the composite system A^ B that 
cannot be constructed in this way is said to be entan- 



gled. Here we will be interested not only in whether or 
not states are entangled, but rather in the amount of 
entanglement in quantum states. Several possible def- 
initions for the amount of entanglement has been sug- 
gested: The entanglement of formation^, the asymp- 
totic entanglement of formation]^ and the asymptotic 
entanglement of distillation [^ |22 1 . All these definitions 
are equally suitable for the purposes of this paper, and 
the results hold for any measure of entanglement which 
satisfies certain continuity requirements. 

To study the behavior of entanglement in noisy 
quantum computers, we define the notion of entan- 
glement length. Roughly speaking, the entanglement 
length is the rate of decay of the entanglement between 
two disjoint sets of particles, as a function of the dis- 
tance between the two sets. This is the analogous quan- 
tity to correlation length in statistical physics, except 
that here we will be interested in correlations between 
two subsets of the system, rather than in two-point cor- 
relations. We study the behavior of the entanglement 
length in the noisy quantum computer as a function of 
the noise rate. We find that there exists a noise rate, 
rji , which depends on the geometry of the system, such 
that the entanglement length is finite for 771 < i] < 1. 
This means that the entanglement between two sets 
of particles decays exponentially with the distance be- 
tween the two sets for this range of noise rates. On 
the other hand, the entanglement length is shown to 
be infinite in the range < ry < ryo, where 770 is the 
threshold for fault tolerance. This means that the en- 
tanglement between two sets of particles is independent 
of the distance between the two sets. These two facts 
show the existence of a phase transition in entangle- 
ment at a non-trivial noise rate 770 < 7?c < 7?i- The 
system in the sub-critical regime behaves quantumly 
even on the macroscopic scale- two sets of particles, 
within macroscopic distance, can share a lot of entan- 
glement, so there is long range entanglement in the sys- 
tem. In the super-critical phase, where entanglement 
decays exponentially with the distance, the system be- 
haves classically on the macroscopic scale, because two 
macroscopic subsets within macroscopic distance are 
practically non entangled. 

The results here are by no means specific to quan- 
tum computers. In fact, we show that macroscopic 
classical behavior is expected above the critical noise 
rate for any macroscopic quantum system with local 
interactions, and local noise, where we make the ad- 
ditional assumption that there is time separation be- 
tween two interactions in which one particle partici- 
pates. This shows that a phase transition in entangle- 
ment is expected in any such quantum system which 
exhibits long range entanglement in the absence of 



noise. Moreover, our results can be verified experimen- 
tally in any such system, as long as the density matrices 
of subsystems can be measured accurately enough. The 
entanglement length can then be numerically computed 
(this is a difficult computational task, but possible for 
small subsystems) and compared to the finite entangle- 
ment length which is predicted by our analysis. 

The emergence of classical macroscopic behavior in 
large quantum systems has been an intriguing area of 
research for the last several decades. Perhaps the most 
common and acceptable explanation so far is by deco- 
herence, i.e. interactions with the environment which 
cause the system to lose its quantum features. See, 
for example, |24, nfl and references therein. This ex- 
planation, however, predicts a gradual transition from 
quantum to classical behavior. The most interesting 
implication of the results presented in this paper is 
that they suggest an alternative way to explain the 
transition from quantum macroscopic behavior to clas- 
sical macroscopic behavior in certain physical systems, 
which is qualitatively different from the standard grad- 
ual explanation. The origin for the abrupt phase tran- 
sition from quantum to classical predicted by our re- 
sults is that in our model, we combine the decoherence 
process with the assumption that noise, as well as inter- 
actions, are local, where the behavior we are interested 
in is global. 

The first part of the proof involves showing that 
the entanglement decays exponentially with the dis- 
tance when 7] is larger than a certain threshold. To do 
this, we use a method due to Aharonov and Ben-Or[Q 
to present the density matrix by mixtures of clustered 
states, and then we study the behavior of the sizes of 
the clusters, evolving in time, using a mapping of the 
problem to percolation. Known results from percola- 
tion theory lljl imply an exponentially decaying bound 
on the probability for distant sets of particles to be 
connected in the same cluster, and this implies expo- 
nentially small entanglement between the two sets. For 
the second part of the proof, i.e. in order to show that 
the entanglement length is infinite for weak noise, we 
use fault tolerant quantum computation, which enables 
one to create long range entangled states in the noisy 
quantum computer. 

We start by defining the notion of entanglement, 
and entanglement length, and then proceed to prove 
the strong noise case and the weak noise case. We con- 
clude with several open questions regarding possible 
implications to the transition from quantum to classi- 
cal. 



2 Entanglement 



The notion of entanglement is associated with a state 
of a quantum system composed of two systems, A and 
B. The term entanglement refers to the quantum cor- 
relations between A and B, in a state which lives in 
A 1^ B. Remarkably, two parts of a composite quan- 
tum system can exhibit very strong correlations which 
cannot be explained classically, unless one drops a very 
important axiom in physics, namely locality. The re- 
markable phenomena which can be exhibited due to 
entanglement between two quantum systems were first 
discovered by Einstein, Podolski and Rosen pi[|, more 
than 60 years ago, and manifested in Bell's inequal- 
ity more than 30 years ago[ffl pi. However, the elusive 
phenomena of entanglement is still far from being un- 
derstood. 

In this paper we will be interested in the amount of 
entanglement in quantum systems; we therefore need a 
good measure of entanglement. One very important re- 
quirement on such a measure is that the entanglement 
in any state cannot increase by classical communication 
and local quantum operations on the A and B sides 
separately; this is the whole essence of the term entan- 
glement. We will denote such a process involving local 
operations and classical communication by LOCC. A 
good measure of entanglement should also be additive. 

A natural way to construct a measure of entangle- 
ment is to ask whether there is an elementary unit 
of entanglement, so that any state can be constructed 
with LOCC given sufficiently many such entanglement 
units. It turns out that there exists exactly such a unit: 
the Beh state, ^(|0)«)|0) + |1)«)|1)). It was shown [H?) 
that any bipartite quantum state can be generated by 
Alice and Bob using only LOCC operations given that 
sufficiently many Bell states are a priori shared between 
A and B. One can try to use the number of elementary 
units required to construct a state as a good measure of 
the entanglement in this state. It is reasonable to take 
the asymptotic limit of such a process, and to define 
the entanglement in a state as the following limit. Let 
(f> be our quantum state, and let fc„ be the number of 
Bell states required to generate a state, (/>„, and let 0„ 
approach the state 0^" as n tends to infinity. Bennett 
et aM defined the infimum over all such processes, of 
kn/n, as n tends to infinity, as the asymptotic entan- 
glement of formation of the state (p. Let us denote this 
measure by _E?°. This measure is clearly additive, and 
can also be shown not to increase by LOCC0. 

An equally natural definition would be the converse 
one, called the asymptotic entanglement of distillation, 
in which one is interested in generating as many Bell 



states as possible by applying LOCC on many copies 
of the original state (p. The asymptotic limit of the 
ratio between the number of Bell states generated in 
this way, and n, was defined in H to be the asymptotic 
entanglement of distillation. A more rigorous definition 
was given by Rains |Q. Let us denote this measure by 

Fortunately, for pure states these two measures co- 
incide, and have a very beautiful form. As was shown 
in |p|, they are exactly the von-Neumann entropy of 
the reduced density matrix on one part of the system. 



E{A:B,\^){cI>\)^S{\cI>){^\\a). 



(1) 



The entropy of entanglement thus possesses both addi- 
tivity and monotonicity under LOCC, and also behaves 
nicely in many other ways. 

The situation for mixed states, however, is much 
more interesting. It turns out that though the asymp- 
totic distillable entanglement and the asymptotic en- 
tanglement of formation coincide on pure states, there 
are very interesting differences between them when 
mixed states are considered. Clearly, the asymp- 
totic entanglement of distillation is not larger than the 
asymptotic entanglement of formation ||7|] . The question 
of whether there exist states in which E?' is strictly 
larger than E°^ is still open. This irreversible pro- 
cess, in which not all of the entanglement which was 
inserted into the state can be distilled, is called bound 
entanglement^^ and is now being extensively studied. 

The asymptotic entanglement of formation is be- 
lieved to be equal to the following quantity, called 
the entanglement of formation^, and denoted by Ef. 
Ef (p) is the least expected entropy of entanglement of 
any ensemble of pure states realizing p, or more for- 
mally: 

Ef{A:B,p)= min Y.w,E{A : B,\a,){a,\). 

^.M>i|Qi)(ai|=p . 

(2) 
The question of whether Ef is equal or not to E'?' 
depends on whether Ef is additive. It is believed that 
indeed it is the case that Ef is additive, but this is not 
known. 

Let us survey what is known about the above three 
entanglement measures, in terms of convexity and con- 
tinuity. Entanglement of formation is trivially con- 
vex. Asymptotic entanglement of formation can also 
be shown to be convex, using the law of large numbers. 
Currently it is not known whether asymptotic distill- 
able entanglement is convex or not. As for continuity, 
the situation is even less clear. It is known that the 
above three entanglement measures are continuous, in 



the sense that if a sequence of density matrices Un con- 
verges to a density matrix p in the trace metric, then 
the entanglement in an converges to the entanglement 
in p: 



lim 



o-ji = P 



lim E{an) = E{p). (3) 



However, we will be interested in how different can the 
entanglement of two very close density matrices be. 
Entanglement of formation was recently shown [gO| to 
have very strong continuity properties, in the following 
sense: Given two density matrices of a bipartite Hilbert 
space of dimension d x d' , which are within e distance 
one from another, in the trace metric, the entanglement 
of formation of the two matrices is at most e times 
some linear function in log(d) and log(d'), plus a term 
independent of d and d' which goes to as e goes to 0: 

\Ef{A:B\p)-Ef{A:B\a)\< (4) 

9|p - a\ log(max{d, d'}) ~ \p ~ a\ log(|p - a\)). 

This strong continuity implies that when two density 
matrices of n finite state particles are polynomially 
close one to another (in the number of particles), the 
entanglement of formation between them is also poly- 
nomially small. It is not yet known whether the asymp- 
totic measures of entanglement possess these nice con- 
tinuity properties, or not. 

In this paper we work with the entanglement of 
formation, Ef, which is known to be both convex and 
strongly continuous. However, it should be stressed 
that the phenomena presented in this paper depend 
very weakly on the exact properties of the measure 
of entanglement which is being used. The results in 
this paper hold, with straight forward modifications, 
for any measure of entanglement E which is contin- 
uous in a sufficiently strong sense, meaning that two 
density matrices which are e apart have entanglements 
not different by more than e times some polynomial in 
the number of particles. 



3 The Model of the Quantum 
System 

We are interested in quantum systems composed of n 
two-state particles, embedded on a d-dimensional lat- 
tice. Such quantum particles are usually called qubits in 
the context of quantum computers. The Hilbert space 
of n such particles is the tensor product of n two dimen- 
sional complex vector space, C^, where the basis of C^ 
is standardly taken to be |0) and |1). The system is ini- 
tialized with a certain state (usually a tensor product 



state, but not necessarily) and evolves in time via inter- 
actions between the particles. Time is discretized into 
time steps and all interactions are assumed to be in- 
stantaneous and occur at integer times. In this model, 
particles interact only with their nearest neighbors on 
the lattice. An important assumption is that one par- 
ticle cannot interact with more than one other particle 
at a time. For simplicity, we will assume that the par- 
ticles interact alternately with particles to each of their 
sides. For one dimension, i.e. an array of particles, this 
means that a particle interacts with a particle to its left 
and to its right alternately. The interaction graph can 
be easily viewed in a d -I- 1 dimensional scheme, which 
for d = 1 looks as follows: 
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Figure 1: The vertical axis corresponds to time, and 
the horizontal axis corresponds to space. Horizontal 
edges connect two interacting particles. Particles in- 
teract alternately with particles to their left or to their 
right. 



Two particles can interact via an arbitrary interac- 
tion, and we do not assume anything about the nature 
or strength of each interaction. After the interactions 
are turned off, and before the next step of interactions 
is turned on, we apply the noise step. The noise is as- 
sumed here to be local and stochastic, meaning that 
each particle with a certain probability -q undergoes an 
arbitrary fault process. As was shown in |p, 13 such d- 
dimensional noisy quantum circuits are capable of per- 
forming fault tolerant quantum computation, as long 
as the noise rate is smaller than a certain threshold. 
The threshold, however, is worse than the threshold 
without the nearest neighbor restriction, by one or two 
orders of magnitude, depending on the dimension. 



We make here another assumption, and restrict the 
noise to be one of the following two processes. The 
first process, namely independent stochastic collapses, 
is a process in which at each time step, each particle 
is measured with independent probability 77, in a fixed 
but arbitrary basis. Alternatively, we can use the de- 
polarization model, in which at each time step, each 



particle, with independent probability rj, is replaced 
by a particle in a completely mixed state. In the rest 
of the paper, we will assume that the noise model is 
independent stochastic collapses, but all results can be 
easily stated using stochastic depolarization. 

It should be noted that the results of the paper 
hold when relaxing several of the assumptions we have 
made. The results apply to particles with any finite 
number of states, not necessarily qubits, and different 
particles need not have the same number of possible 
states. The exact form of the alternating interactions 
is not important, since any interaction graph of nearest 
neighbor interactions is a subgraph of the alternating 
graph, if time is scaled by a factor of 2d. The assump- 
tion of instantaneous interactions is also not essential, 
as long as the interactions last for less than a time 
step, so that there is a time interval in which dccoher- 
ence takes place when the particle is not participating 
in any interaction. The results hold also in the case of 
noisy interactions, with some noise rate 6. We will see 
that the proof of the upper bound on the entanglement 
length, for large i], holds regardless of the amount of 
noise in the interactions, S, because the proof only uses 
the noise occurring between interactions. The proof of 
the lower bound on the entanglement length, for small 
T], goes through as long as rj + S is smaller than the 
threshold for fault tolerance. Hence, it is straight for- 
ward to include in our model the above generalizations. 
For simplicity, however, we will work with the model 
defined above, of two state particles with noiseless in- 
stantaneous interactions. 



4 Entanglement Length 

For quantum systems which are embedded on a lat- 
tice, the notion of distance between sets of particles is 
well defined. In this case, one can define the entangle- 
ment length in the system. We would like to define an 
analogous quantity to the standard correlation length 
from statistical physics. In this case, one says that the 
correlation length in a physical system is ^ if the cor- 
relation between the outcomes of a certain observable 
O measured at sites a and b, decays exponentially with 
the distance between them, d(a, b) , where the distance 
is scaled in the decay factor by ^: 



< OaOt >oc 



(5) 



More precisely, ^ is defined to be the following quantity: 

-log < OaOb > 



In analogy with correlation length, we could have 
defined the entanglement length in the quantum system 
to be /i if the entanglement between two particles, a 
and 6, decays exponentially with the distance between 
them, d(a, 6), where the distance is scaled in the decay 
factor by /i: 

d(a,b) 

E{a:b)(xe ~. (7) 

However, there are a few problems with this defi- 
nition, which will force us to modify it slightly. The 
first modification is necessary due to the fact that en- 
tanglement is a non-local quantity. It might well be 
that the system contains a lot of entanglement, but 
small subsets of the system are completely unentan- 
gled. For example, in fault tolerant quantum comput- 
ers, the entanglement is bound to be shared by large 
sets of qubits, and in order to see entanglement it is 
necessary to probe large subsets of the system. We 
will therefore be interested not in two point correla- 
tions, but in entanglement between two sets A and B 
of arbitrary sizes. 

Another problem is the following. In systems which 
are homogeneous in space and time, one can easily 
take the limit of the size of the system to infinity, and 
therefore the asymptotic behavior in equation (|6|) is 
well defined. However, we are interested in fault toler- 
ant quantum computers, which are not homogeneous 
in space, nor in time. Roughly speaking, we will say 
that the entanglement length in the system is /i if the 
entanglement between any two sets A and B is bounded 
by a function which decays exponentially with the dis- 
tance between the sets, where the decay factor is scaled 
by /x. The fact that we are interested in a bound, and 
not in exact behavior of any pair of sets, allows for non 
homogeneity in space. To allow for non homogeneity 
in time, we will consider the average entanglement be- 
tween A and B over the time from t — to t = oo. 
This corresponds to the following behavior: 



<E{A:B)>Zo^poly{\A\,\B\)e- 



(8) 



C 



lim 

d{a,b)i >oo 



(i(a, b) 



(6) 



where |A| is the number of particles in A and similarly 
for B. We allow the additional polynomial factor due 
to the fact that for sets which are not too large com- 
pared to the distance between them, the exponential 
decay dominates the polynomial in the sizes of A and 
B, and what we will see is merely an exponential de- 
caying behavior. We claim that it is not reasonable to 
consider two sets of particles which are very large, and 
to study the behavior of the entanglement they share as 
a function of the distance between them, in the range 
where that distance is extremely small compared to the 
sizes of the sets. 



The characterization of jj, by equation (|g) is very 
helpful to keep in mind. We can also make the def- 
inition of entanglement length more rigorous, by giv- 
ing it a similar form to that of equation (o). This 
would be useful when one actually wants to calculate 
the entanglement length. In order to do that, we 
first need to make the notion of a quantum system 
more precise. In the non-homogeneous case, it is not 
clear what the notion of an infinite system means. We 
therefore define a quantum (infinite) system to be a 
sequence of quantum systems, Q„, where Qn consists 
of n particles. We think of n as growing to cxi, but for 
a given n Qn is a finite system in space, which evolves 
in time from t = to t = oo. Since each Qn is fi- 
nite in space, in order to take a limit similar to that of 
equation (g), we need to consider a sequence of pairs 
of sets, A and B which belong to larger and larger sys- 
tems. We thus add a subscript n to the subsets An 
and Bn, indicating that they belong to the quantum 
system Qn- We would now like to translate the fact 
that we are interested in sets which are not too large 
compared with their distance to a precise restriction on 
the sequences of sets {An},{B„}. The weakest condi- 
tion which we can impose, to avoid pathologic cases, is 
that lim„, — >oo \An\ ■ \Bn\/exp{d{An, Bn)) = 0, meaning 
that the sizes of A„ and i?„ are not growing exponen- 
tially or faster than exponentially with the distance 
between them. Finally, we want to take care of the 
fact that we are interested in the largest entanglement 
length which can be observed in the system. This cor- 
responds to taking the infimum over all such sequences 
of An and Bn- All this translates to the following def- 
inition: 



Definition 1 The entanglement length fi of a quantum 
system {Qnl^i is defined by: 



u ~ inf liminf 



711^ *00 



-log < EjAn ■- Bn) >t 

d{An,Bn) 



where for all n, An and Bn are disjoint sets in Qn, and 
the sequences {An}'^:^i, {Bnj^^i satisfy liijin: — ►oo \An\- 

\Bn\/eXpidiAn,Bn))^0- 

Note that if one plugs into definition (nl) the exponen- 
tial behavior of equation (0), the contribution of the 
polynomial factor in equation M) tends to zero due 
to the requirements on An and i3„, and the correct /i 
pops out. Though definition ([l|) might seem compli- 
cated, calculating the above infimum turns out to be 
very simple in all our applications. 



5 Clustered Density Matrices 

We now proceed to study the entanglement length in d- 
dimensional noisy quantum circuits, in the strong noise 
regime. In this case, we will try to bound the entan- 
glement in the system from above. Now a very useful 
observation is in place. After a particle was hit by the 
noise process, it is no longer entangled with the rest 
of the system. In other words, in both noise processes 
which we consider, the density matrix after applying 
the noise process can be written in the form 



£p= {l-v)p + V^PtP? ®p1 



(9) 



where the index q refers to the noisy particle, and Q 
refers to the rest of the system. For example, for the 
stochastic noise process, in which the last qubit is mea- 
sured in the basis {|0), |1)}, the resulting density ma- 
trix would be of the form: 

1 
fp=(l-7;)p + 77^Pr(j)pf ®|z)(i| (10) 



i=0 



.Q 



where p^ is the density matrix of all but the last qubit, 
under the condition that the last qubit is measured to 
be in the state |i). We use this observation as follows. 
We will aim to present the density matrix in such a way 
that lack of entanglement translates to tensor product 
structure. In other words, we will present a density 
matrix as a mixture of tensor product states, as follows: 



i 
P^{t) - pl{t)<»---pT'{t) 



(11) 



where p^ (t) is a density matrix which describes a set of 
particles Aj , and for each i the sets Al are a partition 
of the system. These sets of supposedly entangled par- 
ticles are called clusters. It should be understood here 
that given a density matrix, there is no single way to 
present it as a mixture of clustered states. However, 
we will define the representation according to the dy- 
namics of the process which generated p, so that our 
representation will be well defined. 

Our goal would be to find a way to represent the 
matrix as a mixture of clustered states with as small 
clusters as possible. The intuition is that we want to 
give an upper bound on the amount of entanglement 
in the system. When all the clusters are of size one, 
there is no entanglement in the system. We will see 
later, that this can be generalized to say that small lo- 
calized clusters imply no entanglement between distant 
sets. We will thus try to keep the clusters as small as 



possible. The way we do this is as follows. In a quan- 
tum computer, the initial state is a basic state, which 
is a pure state in which all qubits are in tensor product 
with one another: 



p = pil) (g) p{2)(g)---(g) p{n). 



(12) 



Thus, for i = 0, all clusters are of size 1. Given any 
clustered states description of the density matrix at 
time t, we can obtain a clustered state description for 
the matrix at time i + 1 as follows. From each partici- 
pant in the mixture, Pi(t), we obtain pi{t+l) which will 
be a mixture of clustered states. p{t+l) will then be a 
mixture of all pi{t + 1). To obtain pi{t + 1) from pi{t), 
we first apply the interaction step, and then apply the 
noise step. 

To apply the interaction step of time t, we apply for 
each interaction at that time step the unitary transfor- 
mation corresponding to the interaction, on the appro- 
priate pair of particles. If the two particles are from 
one cluster, then we simply apply the appropriate uni- 
tary matrix, corresponding to the interaction, on the 
density matrix describing this cluster, and there is no 
need to change the clusters. However, if the two par- 
ticles are from two different clusters, we can no longer 
keep the two clusters in tensor product, because in gen- 
eral they will be entangled. Therefore, we first join the 
two clusters together, by taking the tensor product of 
the density matrices describing the two clusters, and 
then apply the appropriate unitary matrix on the new 
big cluster. The resulting state after all interactions 
of time step t were applied is therefore, in general, a 
clustered state with larger clusters than the state pi (t) . 

We then apply the noise step on the resulting clus- 
tered state. Recall that a measurement detaches a par- 
ticle from its cluster, and thus after a measurement the 
particle is a cluster of its own. To apply the noise pro- 
cess, we transform the state to a mixture of states, 
which arc the results of all possible combinations of 
which particles where measured, with the appropriate 
probabilities. Clusters in the state can only shrink due 
to this process. 

We would now like to understand the typical size 
of clusters in this representation of the density ma- 
trix. Before we do that in a more formal way, let us 
gain some intuition. If the system were noise free, very 
soon all the clusters would become one giant cluster of 
n particles. What makes the situation more interesting 
are the stochastic collapses, which separate a measured 
particle from its cluster, thereby decreasing the size 
of the cluster by one, and creating another cluster of 
size one. One can view the noisy quantum evolution in 
time as a struggle between two forces: The interactions, 
which tend to join clusters and entangle the different 



parts of the system, and the stochastic collapses, which 
tend to detach particles from their clusters, thereby de- 
stroying this entanglement constantly. A crucial point 
here is that the two competing forces are matched in 
power, since they both operate on a linear number of 
particles 9{n) each time step. We thus expect a crit- 
ical error rate, rjc, at which the two forces are equal, 
and at which some transition between the dominance 
of the entangling interaction process transforms to the 
dominance of the disentangling noise process. We now 
go on to sec this phenomenon more rigorously, using a 
map to a percolation process. 



6 The Percolation Process 

It turns out that the dynamics of the clusters in the 
above description are intimately connected with a per- 
colation process on the quantum circuit itself. The per- 
colation process on the graph is defined as follows: For 
each time step, each vertical edge, along the i'th wire, 
between time t and i + 1 , is erased with independent 
probability rj. In the cluster picture, this corresponds 
to the collapse of the i'th particle between time steps 
t and t + 1, which disentangles the particle's past and 
future. Thus an event in the probability space of the 
noise process, i.e. a specific choice of which particles 
collapsed during the process, is mapped to an event 
in the percolation process, in which the corresponding 
vertical edges are erased. Since events in the stochas- 
tic noise process correspond to members in the mixed 
density matrix, we have a map between clustered states 
arising from our cluster dynamics, and realizations of 
the percolation process. This map preserves the prob- 
ability measure. 

We now claim that clusters in the clustered state 
correspond to connected components in the percolation 
process: 

Lemma 1 Correspondence Lemma: Two parti- 
cles a and b are in the same cluster at time t, in one 
realization of the noise process in the cluster model, 
iff (a, t) and (6, t) are connected in the corresponding 
realization of the percolation model. 

Proof: To prove this combinatorial lemma we use 
induction on t. For the base of the induction, t = 0, 
the correspondence is true by definition. Let us now as- 
sume that the lemma is correct for i, and prove for t+1. 
To apply the induction step, the following observation 
comes in handy. Each path that connects {a,t+l) and 
(6, t -|- 1) in the percolation process, is actually a con- 
catenation of alternating paths, occurring either after 



time t or at times up to t. We denote the points at 
which the different concatenated paths connect one to 
another by (xi,ii), ..., {x2k,t2k)- It is easy to sec that 
there is always an even number of such points, and 
that ti = t2 = ■■■■— t2k = t. Let us call the parti- 
cles xi, ■■■,X2k the connection particles. We shall also 
denote a = xq, b = X2k+i- A schematic example for 
the one dimensional quantum circuit case is shown in 
figure 2. 
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Figure 2: A path connecting the two particles, a and 
b, at time t + 1 can be represented as a concatenation 
of paths which are restricted alternately to the time 
intervals [0, t] and [t, t+1]. The path (a, t+1) i — > (6, t+ 
1) in the figure is a concatenation of the paths {a,t + 

1) I > {Xi,t) I > {X2,t) I > {X3,t) I > ix4,t) I > 

{b,t+l). 

Let us now prove the first direction: Let a and b be 
two particles connected at time i + 1 in the percolation 
model. We want to show that a and b are in the same 
cluster at time i + 1. To see this, we show that all 
particles a ~ x^^xi, ..^X2k,X2k+i = b are in the same 
cluster at time i + 1. For the pairs of particles X2i+i 
and X2i+2i i-e. pairs in which the first particle has odd 
index, this is true since they are connected by a path 
confined to time steps t or earlier, so by the induction 
assumption, they are in the same cluster at time t. 
Moreover, none of the connection particles collapsed 
between time steps t and t + 1, due to the fact that 
they connect between a path before time t and a path 
after time t (or vice- versa). Therefore X2i+i and 2:2^+2 
are in the same cluster also at time t + 1. (In the 
schematic example, this shows that particles xi and X2 
are in the same cluster at time t+1, and similarly x^ 
and Xi.) Now by definition of the connection particle, 
there is a path after time t connecting X2i and X2i+i, 
which means that X2i and X2i+i interact at time t+1. 
At the edges of the chain, i.e. for z = or i = fc, it 
might be that X2i and X2i+i are the same particle, and 
therefore they are trivially connected. (In the example, 
this corresponds to the interaction at time t+1 between 
a and xi, and between X2 and X3, and to the fact that 



b = X4. ) By the definition of the clusters' evolution in 
time, the fact that particles interact imply that their 
clusters are joined, and therefore X2i and X2i+i are in 
the same cluster at time t+1. Combining this with the 
fact that X2i+i and 0:2^+2 are connected at t+1, this 
implies that all the particles a = xq, xi, .., X2kiX2k+i = 
b are in the same cluster at time t+1, which completes 
one direction of the induction step. 

Let us now prove the other direction of the induc- 
tion step. We want to prove that there is a path con- 
necting (a, t + 1) and (6, t + 1) in the percolation pro- 
cess, assuming that a and b are in the same cluster at 
time t + 1 . This cluster of time step t + 1 , which con- 
tains a and 6, was generated by joining together several 
smaller clusters, which existed after the noise step of 
time step t. It is easy to see that there is a subset of 
those clusters, C_i, Co, ..., Ck-i,Ck, such that a <E C_i 
and b € Ck, and such that each two subsequent clus- 
ters, Ci and Ci+i, were connected at time t + 1 by a 
unitary gate. Note that each cluster C'i, except maybe 
C_i and Ck, consists of at least two particles: one 
particle, denoted by X2i+i, participates in an interac- 
tion with a particle from the preceding cluster in the 
chain, Ci_i, and the other particle, a;2i+2 7 participated 
in an interaction with a particle from the next cluster 
in the chain, Ci+i. To construct a path from a to 5 
at time t+1, we first note that by the induction as- 
sumption, there is a path connecting the two particles 
in the same cluster, X2i+i and a;2i+2, at time t. More- 
over, these particles did not collapse between time step 
t and t+1, because if they did collapse, they would have 
belonged to a cluster consisting of one particle exactly. 
Hence there is a path connecting them at time t + 1 . 
A path between the second particle of Ci and the first 
particle of C^+i exists because by definition they inter- 
act at time t + 1. Paths from a to the first particle xi, 
and from the last particle X2k to b, exist either because 
they interact at time t, or they are simply the same 
particle. This enables us to construct a concatenated 
path from a to 6 at time t + 1. □ 

We can now investigate the sizes of the connected 
components in the percolation model and then trans- 
late our findings to the cluster model. Such percola- 
tion processes are known to exhibit a phase transition 
at a critical point, below which there is a connected 
component of linear size, but above which the typical 
connected components is of logarithmic size. More- 
over, the connected components in the super-critical 
phase are localized, in the sense that the probability 
for two particles to be in the same connected compo- 
nent decays exponentially with the distance between 
these particles. 

Let us first release the restriction to a percolation 



on a square of size n x T and consider the infinite lat- 
tice. We show the existence of a phase transition in 
connectivity for percolation on the infinite lattice, and 
from this we will be able to get information about the 
finite case. One more simplification which is useful is to 
notice that by contracting each edge which corresponds 
to an interaction to one point, we do not change the 
connectivity properties of the process, and the resulting 
percolation process is the usual model of percolation. 
For example, in the two dimensional lattice associated 
with the one dimensional quantum circuit, the inter- 
action edges are exactly the horizontal edges in figure 
1, and so after contracting each of these edges to one 
point, the resulting percolation process is the standard 
percolation on the square lattice, which is rotated by 
45 degrees. The contraction therefore transforms the 
problem to standard bond percolation on translational 
invariant lattices, which is pretty well understood. 

In bond percolation on translational invariant lat- 
tices, one usually uses p as the probability for one edge 
to be present, so in our case p = I — r]. One can de- 
fine the critical p, pc, to be the smallest probability 
in which the point belongs to an infinite connected 
component with positive probability. In translational 
invariant lattices, this is also equal to the smallest p in 
which the expected size of O's connected component is 
infinitc|l|, 0, 1 

p, = inf{p I Prp{\H{0)\ = (^) > 0} = (13) 
= inf{p\Ep{\HiO)\)^^}. 

where H{Q) is the connected component of 0. 

A theorem by Hammersleyg, |lj], asserts that for 
translational invariant lattices, for p < Pc, the prob- 
ability T{x,y : p) for x to be connected to y, decays 
exponentially with the distance: 



Lemma 2 i < Pc{d + 1) < ^ro 



t(x, y : p) < exp 



d{x,y) 



(14) 



with ^(p) < CX3. Above Pc, the probability for to be 
connected to infinity is larger than by definition of 
Pc, so ^{p) = oo for p > Pc- 

Let us denote hy pc{d + \) the critical pc for bond 
percolation on the lattice corresponding to a quantum 
circuit of dimension d. For one dimensional quantum 
systems, it is easy to calculate the critical probability 
Pc(l + !)• After contracting the interaction edges, we 
simply get percolation on Z^, for which it is known [p^ 
that the critical probability is a half. Hence: 

Pc{l + l) = \. (15) 

For higher dimensions, we can bound Pc{d + I) away 
from and 1, but the bounds are not tight: 



Proof: The upper bound comes from the fact that the 
projection of the d-|- 1 percolation with parameter p on 
a 1 -|- 1 percolation gives percolation in I -I- I with pa- 
rameter p"*. This is true since after a particle interacts 
with a particle along one axis, it waits d time steps be- 
fore it interacts again with a neighbor along the same 
axis. This gives the upper bound, since if the origi- 
nal process had exponentially decaying correlations, it 
cannot be that in the projected process we are above 
the phase transition where is connected with constant 
probability to infinity. The lower bound is derived from 
a standard argument which reduces the problem to a 
branching process. A branching process is a process 
which starts with one node, and each nodes gives birth 
to k nodes with some probability distribution p(fc), in- 
dependent of the other nodes. It is a standard result 
(See [E2I , for example) that when the expected number 
of descendents for each node is less than I, the dynasty 
dies in finite time with probability 1. To construct 
the corresponding branching process, observe that the 
degree of the interaction graph, after contracting the 
horizontal edges, is exactly 4, regardless of d. Starting 
with the point 0, we regard each of its neighbors to 
which it is connected in the percolation as an ancestor 
of a dynasty. The descendents of each such node are 
all its neighbors, except for 0, in the percolation pro- 
cess; Each descendent has its own descendents, and so 
on. When we encounter a node which is already in the 
dynasty, we do not count it again; this way we have 
a tree. Clearly, if the branching process is finite, then 
the connected component of in the percolation pro- 
cess is definitely finite. However, each neighbor of is 
the ancestor of a dynasty in which the expected num- 
ber of descendents of each node is exactly the number 
of its neighbors which are not yet in the dynasty, times 
p. Since the degree of the graph is 4, and one of the 
neighbors is the node's ancestor, the expected number 
of descendents is at most 3p, which is less than 1 for p 
less than i. This gives the desired result. □ 

The analysis of the percolation process has taught 
us that when p is smaller than the critical point for 
percolation, pc, the connected components in the sys- 
tem are small. Going back to the density matrices, 
using the correspondence lemma, this implies that for 
7] < 1 — Pc, the clusters are small and localized. This 
will be used in the next section to prove an upper on the 
entanglement between distant sets in this noise regime. 



7 Finite Entanglement Length 

From the correspondence to percolation, we have that 
for r/ < 1 — pc, the density matrix of the quantum sys- 
tem can be approximated by a mixture of clustered 
states with localized clusters of logarithmic size. Thus, 
distant subsets of particles are with high probability 
contained in non-intersecting clusters, i.e. most of the 
weight of the density matrix is concentrated on states 
in which there is no entanglement between the two sub- 
sets. The weight of the states in which there is entan- 
glement between the two sets decays exponentially with 
the distance between the sets. By continuity of entan- 
glement, this implies that the entanglement between 
the two sets decays exponentially with the distance. 

We can now show, that the entanglement between 
any two sets of particles decays exponentially with the 
distance between them, when the noise rate rj is such 
that 1 — 77 is sub-critical in the percolation process. The 
rate of the decay is the entanglement length of the sys- 
tem. The entanglement between the two sets becomes 
negligible already when the distance is of the order of 
log(n) particles. This translates to the following theo- 
rem; 

Theorem 1 Consider a d dimensional quantum, cir- 
cuit with nearest neighbor interactions, subjected to lo- 
cal noise of the type of stochastic depolarization or 
stochastic collapses, with noise rate rj. If the circuit 
is initialized with an unentangled state, i.e. a tensor 
product state, and if rj > 1 — pc{d -\- 1) , then the entan- 
glement of formation between any two sets of qubits A 
and B at any time t > decays exponentially with the 
distance between the two sets: 



Ef{A: B) < m\n{\A\,\B\}\A\ ■ \B\e W~^ . 

For a general initial state, a similar formula is true 
except for a correction term which decays exponentially 
with time: 



'fiA : B,pit)) < min{|A|, \B\} {\A\ ■ \B\e~^^- 
+nmm{\A\, |-B|}e"«<i-") 



E 



Proof: Let us start with the simple case, in which 
the initial state is a complete tensor product, i.e. all 
clusters are of size 1. By equation |lj, the probabil- 
ity for two particles A and B to be connected decays 
exponentially in the distance between them. The cor- 
respondence lemma (|l|) implies that the probability 
for two particles to be in the same cluster at time 
t is equal to the probability they are connected in 



the percolation model at time t. i.e. the probabil- 
ity for two particles from A and B to be in the same 
cluster is bounded above by exp{—d{A,B)/(,{l — ry)). 
Thus, the probability for any pair of particles from A 
and B to be in the same cluster is bounded above by 
\A\-\B\exp{-d{A, B)/^{l-ri)). The density matrix can 
thus be written as a mixture of one density matrix with 
weight smaller than \A\ ■ \B\exp{-d{A,B)/^{l - 77)), 
and another density matrix which is a mixture of den- 
sity matrices, where in all these matrices all the parti- 
cles in A are in different clusters than all the particles 
in B. The reduced density matrix to A, B of the second 
matrix is thus separable, and contains no entanglement 
between A and B. By convexity of entanglement of for- 
mation, the entanglement in the entire density matrix 
is bounded above by the entanglement in the first den- 
sity matrix, times the weight of this matrix. The en- 
tanglement of the first matrix is at most the number of 
qubits in the system, and this gives the desired result. 
(For measures of entanglement which are not convex, 
but strongly continuous, one should replace the term 
min{|A|, \B\} by the appropriate polynomial from the 
continuity bound.) 

We now proceed to the general initial state. We will 
give an upper bound for the case in which the initial 
state is one big cluster, and any other case is trivially 
implied by it. To do this, we have to understand where 
we have used the fact that the initial state is not en- 
tangled. This was used for the base of the induction 
in the correspondence lemma, where the fact that all 
clusters are of one qubit corresponds to the fact that 
in the percolation graph, the initial connected compo- 
nents at time are all of size 1. To adapt the situation 
to the case in which all particles are in one big cluster 
at time 0, we add a horizontal line of length n con- 
necting all particles to one big connected component 
at time i = 0. The correspondence lemma then goes 
through. However, equation nj no longer holds. To 
correct it, we add to it a term which corresponds to 
the probability for A to be connected to B by a path 
that goes through time i = 0, i.e. through the addi- 
tional new line we have added to the graph. For such 
a path to exist, both A and B need to be connected to 
time 0. The probability for any one of the qubits in A 
to be connected to any one of the n qubits at time 
is at most n\A\ times the probability for one qubit at 
time t to be connected to one qubit at time 0, which is 
at most exp{— f^l., ■) ) by equation ^ The same argu- 
ment applies for the connection from B to time 0, and 
this gives the desired result. □ 

This shows that the system cannot create entangle- 
ment between far sets of particles: Roughly speaking, 
the typical range of entanglement is microscopic. This 
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is true for any initial condition, where the relaxation 
time to the typical unentangled state is of the order of 
log(n) steps. 

This result implies an upper bound on the entangle- 
ment length in the quantum system above the critical 
noise rate, and in particular shows that it is finite. This 
is done by simply taking the limit in the definition of 
entanglement length (m, which gives: 

Corollary 1 The entanglement length ^{rf) of a d di- 
mensional quantum circuit with nearest neighbor inter- 
actions, subjected to local noise of the type of stochastic 
depolarization or stochastic collapses, with noise rate rj, 
satisfies: 

and in particular fi{ri) is finite for rj > 1 ~ pdd + 1). 

This gives a bound on the entanglement length, in 
terms of the correlation length in classical bond perco- 
lation. The correlation length of a given lattice can be 
easily estimated by computer experiments, and analyt- 
ical bounds are given in m, nj] . 



8 Infinite Entanglement Length 

We now want to concentrate on the other noise regime, 
and show that below the critical noise the entangle- 
ment length is infinite. One might naively think that 
this can be deduced from the fact that the density ma- 
trix is a mixture of clustered states with linear sized 
clusters. However, there is a difficulty in pursuing the 
connection between clusters and entanglement for this 
purpose, because of the following reason. The density 
matrix is actually a mixture of many clustered states. 
The mixture of two clustered states, with very large 
clusters, can be a density matrix in which the clusters 
are of size one. One example of such a case is a mixture 
of the two states, ^(|0") + |1")) and ^(|0") - |1")), 
the mixture of which is a non entangled state. Thus, 
the sizes of the clusters can be used for upper-bounds 
on entanglement, but it is not clear how to use them 
in order to show a lower bound on the entanglement in 
the system. 

We therefore need to use different techniques for 
lower bounds on entanglement. We will use techniques 
from quantum computation. A quantum computer em- 
bedded on a lattice is a special case of the quantum 
systems wc are discussing. The particles are quantum 
bits, and the interactions are fixed according to the 
algorithm. Therefore, corollary |] shows that the en- 
tanglement length is finite above the critical noise rate 



also in fault tolerant quantum computers. For fault 
tolerant quantum computers we can also analyze the 
other side of the noise scale, and show that the entan- 
glement length in the system is infinite if the noise rate 
is below a certain threshold. We will use the threshold 
result m |l^, |l^ for fault tolerant quantum computa- 
tion, which shows that quantum computation can be 
made robust to noise, using quantum error correcting 
codes, as long as the noise is smaller than a certain 
threshold. In fact, here we need the slightly stronger 
version of the threshold result 0, ITSJ, which asserts that 
this can be done even when the quantum system is em- 
bedded on a d dimensional lattice. The threshold is 
then r]o{d), which for d = 1 is estimated to be 10^^ B. 
In the fault tolerant range, two distant sets of qubits 
can be entangled, and remain entangled for a long time, 
with the amount of entanglement independent of the 
distance. 

We now give an example of a quantum computer 
which exhibits entanglement among far parts of the 
system when 77 < 770 (rf), but the entanglement length 
is finite for noise 77 > 1 — Pc{d + 1). The idea is that 
a fault tolerant computer can simulate any quantum 
state, including states which contain entanglement be- 
tween sets of qubits which are far apart. Hence, we 
will construct a quantum algorithm in which there is 
entanglement between two far parts of the system, and 
make it fault tolerant. This can be done in many ways, 
but here is a simple example, for d = I. Divide the set 
of qubits to three sets. A, B, C. We will create entan- 
glement among A and B, while leaving the qubits in 
the middle, C, in a basic state. This will be done by 
constructing the state: 



^(10- 



)a ® |0")c ® |0«)i3 + |1")a ® |0")c ® \1')b) 

(16) 

on a fault tolerant quantum computer, and keeping this 
state for a long time, by applying error corrections. 
This state indeed contains entanglement between the 
two registers A and B, which are n sites apart. The 
algorithm which constructs such a state is very simple, 
and uses only two basic quantum gates: The Hadamard 
gate, which is a one qubit gate applying the following 
unitary transformation 



|0) 



ID 



71<'"' 



|i» 



(17) 



;^(|0)-|1», 



and the controlled NOT gate which is a two qubit gate 
applying the following unitary transformation (said to 
be applied from the first qubit to the second one): 



|a)(g) \b) I — > \a) (g) |a©6), 



(18) 



11 



where © means addition mod 2. Using these gates, it 
is easy to create the state 



1 



— (|0"+^) + |l™+')) 



(19) 



on the first m + q qubits, by applying a Hadamard 
gate on the first qubit and then controlled NOT gates 
from the first qubit to the second, from the second 
to the third, and so on. Then, we want to swap the 
m + 1, ..., m + q qubits to register C. To do this, we 
first swap the last qubit in B with qubits to its right 
until it gets to the last site in C; In the same way we 
bring the one before last qubit in B to the one before 
last site in C, and so on until all qubits in B are in the 
right most sites of C, which achieves the desired state 
with only nearest neighbor interactions. 

This algorithm by itself is not fault tolerant, and 
in the presence of any amount of noise, i.e. rj > 0, 
the entanglement in the system will be lost immedi- 
ately. However, we can make this algorithm fault tol- 
erant by the methods in ||, |l^, [T7|, as long as 77 is 
smaller than rj{d), the threshold for fault tolerance for 
d-dimensional quantum computers Isl O]. These re- 
sults are too complicated to explain here in details. 
In a nutshell, fault tolerance is achieved by encoding 
the qubits using quantum error correcting codes, and 
computing the algorithm on the encoded states, while 
applying quantum error correction on the state fre- 
quently. Each qubit is replaced by polylog(n) qubits, 
encoding its state. The state |0) is encoded by the state 
l^o) of poly log (n) qubits, and similarly |1) is encoded 
by the state 15*1). Let us denote by A' , B' and C" the 
qubits encoding the original sets of qubits A, B, and 
C, respectively. If no fault occurs, at the end of the 
algorithm the state of the system will be in the state 
( |l6|) encoded by the quantum error correcting code: 



I^™)a' ® \SS)c" ® l^o') 



I ST 



E> \SS)c' ® \S1)b' 
(20) 
normalized by a factor of -j^. The entanglement in 
this state will remain there for ever if errors do not oc- 
cur. However, errors do occur. By fault tolerance, this 
means that at the end of the computation the den- 
sity matrix is polynomially close to a density matrix 
p which can be corrected to the correct state (EQ) by 
noiseless quantum error corrections. Due to continuity 
of entanglement, it suffices to argue that such p con- 
tains a constant amount of entanglement. But this is 
true since we know that p can be corrected to the state 
( P0|) by local operations not involving interactions be- 
tween A' and B' . Since entanglement cannot increase 
by local operations, the entanglement between A' and 
B' in p is at least as that in the state encoding (po|), i.e. 
one entanglement unit. The distance between the ac- 



tual density matrix and a correctable density matrix p 
is, by H, at most the number of time steps t divided by 
a polynomial in n. This distance is smaller than some 
constant e as long as the number of time steps is poly- 
nomial in the size of the system n. Thus, by strong 
continuity the entanglement between A' and B' will 
remain bounded from below by a constant for poly- 
nomially many time steps. After polynomially many 
steps, we can replace all qubits by qubits in the state 
|0), and run the whole algorithm again. The average 
entanglement over time from to 00 is very close to 
one, since the time it takes to construct the state is 
much smaller than the polynomial time for which the 
entanglement remains in the system. This proves the 
existence of a non trivial sub-critical side of the phase 
transition: 

Theorem 2 The entanglement length in the d dimen- 
sional fault tolerant quantum computer defined above 
satisfies 

p.{ri) = 00 

for 77 smaller than the threshold for fault tolerance in d 
dimensional quantum computers, i.e. rj < riQ{d). □ 



9 Other Quantum Systems 

The model of a noisy quantum computer actually holds 
not only for quantum systems designed to serve as com- 
putational devices, but for a much broader class of 
physical systems as well. We first claim that putting 
aside the noise process, any d dimensional quantum 
system, in which the particles are located in space with 
low enough density, and in which interactions occur 
only between particles which arc not too far apart, can 
be modeled by a quantum circuit. This can be done 
by discretizing the medium to very small cells, such 
that each cell contains at most one particle. Time will 
be discretized to sufficiently small intervals such that 
a particle can only move to a neighbor cell in one time 
step. Then, the movement of particles can be mod- 
eled by an interaction between an occupied and an un- 
occupied cell, and since the density of particles is low, 
one particle never interacts with more than one other 
particle at the same time, so the notion of quantum 
gate is appropriate. 

The noise model which is used in this paper is quite 
general as well, when low density or instantaneous in- 
teractions systems are considered. During the time in- 
terval in which a particle does not participate in any 
interaction, stochastic collapses are actually equivalent 
to a process of local decoherence. Assume that each 
particle interacts with its own independent thermal 
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bath, and the Markovian assumption is apphed, so that 
the environment of each particle is renewed each time 
step. This corresponds to the process in which the off 
diagonal terms in the density matrix of one particle 
decay by some factor between two interactions [^ . If 
the decoherence process operates for time At, the {i,j) 
element of the density matrix, transforms to 

Pi,j I — > Pi,j exp(-7Ai(l - (5ij)) (21) 

If we set exp(— 7Ai) = 1 — 77, we get 



Pi. 



(1 - i])pi^.j + ripijSi, 



(22) 



which is equivalent to a measurement with probability 
77. Similarly, the depolarization process can be pre- 
sented as a gradual change of the density matrix of 
one particle. 

The above arguments show that the model of noisy 
quantum circuits which we are discussing is interesting 
as a representative of the class of quantum systems 
with macroscopically many finite state particles, local 
instantaneous interactions and local decoherence noise. 
The analysis done in this paper regarding upper bounds 
on the entanglement length in the super-critical phase 
goes through, and therefore theorem n and corollary n^ 
can be generalized to this case. 

In such quantum systems, our analysis provides an 
explanation to the emergence of macroscopic classical 
behavior above the critical noise rate, as will be dis- 
cussed in the conclusions. 



10 Experimental Verification 

Unfortunately, we do not yet have a physical realization 
of a quantum computer of more than several qubits, on 
which the existence of a phase transition in entangle- 
ment length in fault tolerant quantum computer can be 
verified. However, for the more general case of quan- 
tum systems with local interactions and local noise, 
satisfying the requirements of section 0, the bound on 
entanglement length in the super-critical regime can in 
principle be experimentally testable. What is needed 
to perform such an experiment is to be able to mea- 
sure with high enough accuracy the joint density ma- 
trix of two subsets of the system. The entanglement 
between the two sub-systems can then be numerically 
approximated, using a (doable, but extremely difficult) 
minimization over equation g. The entanglement can 
then be found as a function of the distance between 
the sets, from which the entanglement length can be 
deduced. An extremely interesting open problem is 
to give a concrete design for such an experiment, for 



an existing quantum system, and to compare the out- 
comes with the entanglement length predicted by the 
percolation process. 



11 Quantum- Classical 
Transition 

The results presented here suggest that the emergence 
of classical macroscopic phenomena in large quantum 
systems can be attributed, in certain cases, to the fact 
that the noise rate is larger than a certain critical point, 
so that the entanglement length is finite. However, we 
merely introduced in this paper a new phenomena. The 
list of questions which remain open is extremely large, 
and varies on different physical fields. 

The first and most basic question should be how 
general these results are. In this paper, we have been 
able to show the existence of a non trivial sub-critical 
phase in fault tolerant quantum computers. Are there 
more natural quantum systems, in particular systems 
which are homogeneous (or periodic) in space and in 
time, which have local interactions and local noise, 
which are able to maintain long range entanglement in 
the presence of weak or zero local noise for a long time? 
Does a random quantum system, i.e. in which random 
interactions are applied, exhibit long range quantum 
correlations? Such systems will provide more exam- 
ples for quantum systems in which a phase transition 
in entanglement length occurs. 

It should be noted here that the notion of zero noise 
rate does not trivially coincide with that of zero tem- 
perature, i.e. long range entanglement in the ground 
state of the Hamiltonian of the system. An important 
observation is that the model discussed here deals with 
non-equilibrium quantum systems. The quantum sys- 
tems we consider here can be in a steady state, but the 
density matrix is not in the Gibbs distribution of the 
eigenvectors of some Hamiltonian. The reason for this 
is that we did not allow the system to approach equi- 
librium. Our noise model, or the interaction with the 
thermal bath, is limited to local interactions. This is a 
crucial ingredient that causes the phase transition. It 
is the fact that two forces of even power compete: lo- 
cal interactions in the system against local interactions 
with the environment, that gives the critical noise rate. 
The fact that the quantum computer does not achieve 
equilibrium despite the noise is explained by the fact 
that the system is cooled constantly by quantum er- 
ror corrections. It is left as an open problem to further 
investigate possible equilibrium phase transitions in en- 
tanglement, and the connection to the non-equilibrium 
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phase transition presented here. The reader is referred 
to |l^ and references therein for an introduction to 
non-equihbrium phase transitions. We view the fact 
that we did not allow the system to achieve equilib- 
rium as very important in the derivation here. 

Indeed, it is worthwhile to ask which of our assump- 
tions regarding the properties of the physical system 
are essential, and which are technical. Intuitively, the 
locality of both noise and interactions seems crucial. 
This is true since this locality is what makes the two 
competing forces, the interactions which tend to entan- 
gle the system and the noise which tends to disentangle 
it, comparable in power, which gives rise to the phase 
transition. It seems, however, that the assumption on 
the exact form of the noise process might not be so 
important, and neither is the discretization of the in- 
teractions. An important open problem is to relax the 
assumptions used in this paper, and to generalize the 
results presented here to other noise models, and to the 
continuous case. In particular, it is not clear how to 
generalize the results to the case in which the parti- 
cle interacts simultaneously with all its neighbors and 
the environment. It seems that a considerably different 
approach would be needed in this case. 

If the phenomena of the phase transition in entan- 
glement is indeed general, its effect on our understand- 
ing the transition from quantum to classical physics 
needs to be deeply understood. One important ques- 
tion is whether there exists some classical or quasi- 
classical description of the behavior of a quantum sys- 
tem in its super-critical phase. Another, related, ques- 
tion is whether the existence of a phase transition in 
entanglement induces other quantum phase transitions 
at the same critical point in the same system. 

A set of open questions regarding the phase transi- 
tion comes from statistical physics. For example, what 
are the critical exponents related to this phase tran- 
sition? What is the universality class of this phase 
transition? In fact, it is not clear that there is only 
one critical point here. In the case of the quantum 
computer, or other quantum systems, there might well 
be an intermediate regime of noise, in between the two 
thresholds, for which the entanglement behaves in a 
different way, i.e. its dependence on the distance is 
neither an exponential decay nor constant. The ques- 
tion of showing that there is only one critical point at 
which a transition occurs from exponential decay to 
independence on the distance remains open. 

A very interesting problem is to come up with a 
better order parameter related to entanglement, rather 
than the entanglement length. There are many prob- 
lems with the entanglement length as an order param- 



eter. The most important one is that it might be that 
the system is very entangled, but the entanglement be- 
tween two distant subsets is zero. Such is for example a 
system in the state ^jdO") + ll")), for which any sub- 
system is non-entangled. Entanglement in such very 
entangled quantum systems will not be detected when 
looking at sub-systems, and the entanglement length 
will therefore contain no information about the actual 
behavior of the entanglement in the system. Another 
motivation for this question is provided by B, where 
the sizes of the clusters is analyzed relaxing the as- 
sumptions of nearest neighbor interactions. The sizes 
of the clusters in this case indeed transform from log- 
arithmic to linear at a critical noise rate. However, 
the notion of entanglement length cannot be defined 
in a system without geometry, so it is not clear how 
to define an order parameter which exhibits the phase 
transition in this case. 

To summarize, we have discovered a phenomena of 
a phase transition in entanglement in quantum com- 
puters, and in general in quantum systems with lo- 
cal decoherence and local interactions which are able 
to generate long range entanglement in the absence of 
noise. The suggestion to explain the transition from 
quantum to classical macroscopic behavior as a phase 
transition in entanglement is fundamentally different 
from the standard point of view of gradual transition, 
usually explained by decoherence. Our results have ex- 
perimental implications, and raise a long list of open 
problems related to the foundations of quantum me- 
chanics, as well as to quantum statistical physics. 
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